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Рассмотрены алгоритмы определения параметров утечки в трубопроводах, полученные по уравнениям гидродинамических про- 
цессов в них, возникающих при утечках, приведены результаты исследования предложенных алгоритмов определения параме- 
тров утечки. 


Для транспорта нефти и газа актуальной являет- 
ся проблема обеспечения надёжности магистраль- 
ных трубопроводов. При возникновении утечки по 
любой причине необходимо в кратчайшие сроки 
обнаружить сам факт существования утечки и опре- 
делить место её возникновения. Сокращение вре- 
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мени от момента возникновения утечки до момента 
обнаружения с установлением её координаты по- 
зволяет максимально быстро принять меры по ло- 
кализации аварии, сокращая её масштабы и затраты 
на ликвидацию. Для удовлетворения перечислен- 
ным требованиям трубопроводный транспорт дол- 
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жен быть оснащён автоматически работающими 
системами обнаружения утечек. 

В представленной работе рассматривается ме- 
тод определения утечки по изменению профиля 
давления вдоль трубопровода, который базируется 
на фиксации понижения давления в случае возни- 
кновения утечки на контролируемом участке [1,2]. 
Предполагается, что рассматриваемый трубопро- 
вод сплошной, без дополнительных отводов, по ко- 
торому протекает однородная жидкость. 

При разработке алгоритмов использовалась ги- 
дродинамическая модель, описывающая процессы 
утечки в трубопроводе в виде дифференциального 
уравнения и граничных условий [3]: 
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ка; Р - давление в магистрали; 1 - время; х - коор- 
дината по длине трубы; /, сі — длина и внутренний 
диаметр трубопровода; Р - площадь сечения трубы; 
С - скорость распространения волны давления; 
С Т - массовый расход утечки; Я - коэффициент со- 
противления трубопровода; м> - скорость движения 
нефтепродукта; ^ - координата места утечки. 

Для приведённых уравнений может быть полу- 
чено, методом разделения переменных (методом 
Фурье), точное решение [3]: 
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где Р н , давления, соответственно, в начале и на 
конце участка трубопровода. Эти значения давле- 
ний в процессе функционирования трубопровода 
должны поддерживаться постоянными. 



Рис. 1. Схема трубопровода при наличии четырёх датчиков 


Будем рассматривать прямолинейный трубо- 
провод, пролегающий между двумя насосными 
станциями. На каждой станции имеется датчик да- 


вления, а на всем протяжении трубы располагают- 
ся два датчика давления 2), и Д (рис. !)• 

Тогда расчётные формулы для определения па- 
раметров утечки по профилю давления, определя- 
емому по четырём датчикам, будут иметь вид [4]: 
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Р в | и Р т - давление на первом и втором датчиках. 

Как видно из формул (1), решение содержит 
стационарную и нестационарную составляющие. 
Алгоритм (2) основан на стационарном решении 
гидродинамической модели, определяющей про- 
филь давления по четырём датчикам. Проверка 
адекватности модели и работоспособности полу- 
ченного алгоритма была проведена на имитацион- 
ных моделях, разработанных в программном паке- 
те МаіБаЪ, с использованием приложения МаіБаЬ 
- 8 і гті и I і п к . Расчёты проводились для реального 
объекта по обобщённой схеме, рис. 2. 



Рис. 2. Обобщённая схема вычисления параметров утечки 



В результате моделирования получен график 
профиля давления, представленный на (рис. 3). Из 
графика видно, что в месте прорыва трубопровода 
(^=18-10’м) происходит надлом кривой, параметры 
которого могут быть вычислены по формулам (1). 
Это доказывает эффективность алгоритма опреде- 
ления параметров утечки по профилю давления. 


71 







Известия Томского политехнического университета. 2006. Т. 309. № 7 



Рис. 4. Г оафики процесса вычисления координаты утечки и массового расхода при использовании алгоритма по профилю да- 
вления с четырьмя датчиками 


Однако при использовании данного алгоритма 
определение параметров начинается после завер- 
шения переходных процессов в трубопроводе, что 
увеличивает время определения. 

Для повышения быстродействия определения 
параметров утечки с определённой точностью 
можно использовать рассмотренный выше алго- 
ритм обнаружения по профилю давления, опреде- 
ляемого четырьмя датчиками давления, не дожида- 
ясь завершения переходных процессов в трубопро- 
воде при возникновении утечки. 

При этом расчёт выполняется по тем же форму- 
лам, но в модели профиля давления в трубопрово- 
де учитывается нестационарное решение. 

Графики процесса вычисления координаты 
утечки и массового расхода при использовании ал- 
горитма по профилю давления с четырьмя датчика- 
ми с учётом переходных процессов в трубопроводе 
представлены на рис. 4. 

По приведённым графикам видно, что при ис- 
пользовании рассмотренного выше метода на ос- 
нове показаний четырёх датчиков давления для 
определения координаты утечки и расхода с одно- 
процентной погрешностью происходит за время 
5,1 и 4,0 с соответственно, что на порядок меньше 
времени окончания переходных процессов в трубо- 
проводе. 

Для ещё большей оперативности определения 
параметров утечки был разработан алгоритм [5], 
который предполагает применение шести датчи- 
ков, определённым образом расположенных по 
трубопроводу (рис. 5). Тогда выражения для расчё- 
та расхода и координаты утечки имеют вид: 
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Моделирование проводилось для того же объекта. 
Графики зависимости координаты утечки и массового 
расхода при использовании алгоритма по профилю 
давления с шестью датчиками представлены на рис. 6. 
При этом прорыв трубопровода происходит между 
третьим и четвёртым датчиками, то есть <^=18Т0 3 м. 
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Рис. 5. Схема расположения датчиков на трубопроводе при 
шести точечном контроле 


Из приведённых графиков видно, что при ис- 
пользовании метода на основе показаний шести 
датчиков решения для определения координаты 
утечки и расхода становятся точными за время 1 с. 

Были проведены исследования алгоритмов с 
учётом зоны нечувствительности датчиков (Д), а 
именно для /1=5 кПа. Результаты исследований для 
расчёта массового расхода и координаты утечки 
для алгоритма по четырём и шести датчикам при- 
ведены на рис. 7. Данное исследование показало, 
что наличие зоны нечувствительности влияет на 
точность определения параметров утечки, причём, 
большей точностью обладают вычисления при ис- 
пользовании алгоритма с шестью датчиками. 
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Рис. 6. Г оафики процесса вычисления координаты утечки (слева) и массового расхода (справа) при использовании алгорит- 
ма по профилю давления с шестью датчиками 



а Ь 

Рис. 7. Г оафики процесса вычисления координаты утечки и массового расхода при использовании алгоритма по профилю да- 
вления с р=5 кПа: а) при четырех датчиках, Ь) при шести датчиках 


При исследовании представленных выше алго- 
ритмов определения параметров утечки было ис- 
следовано влияние некоторых параметрических 
возмущений, в частности, влияние неточности 
определения скорости звука в жидкости. 

Исследования показали, что при задании в ал- 
горитмах вычисления скорости звука с погрешно- 
стью в 10 % параметры утечки будут вычислены с 
погрешностью в 3 %. Для практического использо- 
вания рассмотренных в работе алгоритмов требует- 
ся более полное исследование влияния на них все- 
го спектра параметрических возмущений. 


Выводы 

1. Получены алгоритмы определения параметров 
утечек по точным решениям уравнений гидро- 
динамических процессов в трубопроводах при 
утечке. Показана в результате имитационного 
моделирования работоспособность этих алго- 
ритмов. 

2. Выбран метод исследования алгоритмов, про- 
ведено моделирование процессов определения 
параметров утечки по трём алгоритмам: по про- 
филю давления без учёта и с учётом переходных 
процессов в трубопроводе, определяемого че- 
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тырьмя датчиками давления и шестью датчика- 
ми давления, определённым образом устано- 
вленными по длине трубопровода. 

3. Сравнение исследуемых алгоритмов показало, 
что при одинаковой точности определения па- 
раметров утечек наибольшей оперативностью 
(быстродействием) обладает алгоритм опреде- 
ления параметров утечки по профилю давления 
с шестью датчиками, установленными опреде- 
лённым образом по длине трубопровода. 

4. В результате исследования влияния нечувстви- 
тельности датчиков давления на результаты 
определения параметров утечки, было выявле- 


СПИСОК ЛИТЕРАТУРЫ 

1. Барбашов Е.Д., Коняхин А.Н. Результаты испытаний параме- 
трической системы обнаружения утечек // Комплексные ре- 
шения в автоматизации для динамического развития нефтега- 
зовой отрасли: Матер. II научно-практич. конф. - 1-2 ноября 
2001 г. - Томск: Изд-во Томского ун-та, 2002. -С. 62-67. 

2. Гольянов А.А. Анализ методов обнаружения утечек на трубо- 
проводах // Транспорт и хранение нефтепродуктов. - 2002. - 
№ 10-11. -С. 5-14. 

3. Степанченко Т.Е., Шкляр В.Н Исследование математической 
модели процесса утечки жидкости в трубопроводе // Совре- 
менные техника и технологии: Труды VII Междунар. научно- 


но влияние величины зоны нечувствительности 
на точность определения. Установлено, что ал- 
горитм с шестью датчиками позволяет вычи- 
слить координаты утечки и её массовый расход 
с большей точностью, чем алгоритм с четырьмя 
датчиками при одинаковом значении зоны не- 
чувствительности датчиков. 

Для практического использования представлен- 
ных алгоритмов, в дальнейшем планируется сделать 
доработку алгоритмов определения параметров 
утечки для участков трубопровода с уклонами и от- 
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большего числа параметрических возмущений. 
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